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Abstract 

We address the perturbative renormalization of massive lattice fermions. We 
derive expressions — valid to all orders in perturbation theory and for all values 
of the bare fermion mass — for the rest mass, the kinetic mass, and the wave- 
function renormalization factor. We obtain the fermion's self energy at the 
one- loop level with a mass-dependent, 0{a) improved action. Numerical re- 
sults for two interesting special cases, the Wilson and Sheikholeslami-Wohlert 
actions, are given. The mass dependence of these results smoothly connects 
the massless and infinite-mass limits, as expected. Combined with Monte 
Carlo calculations our results can be employed to determine the quark masses 
in common renormalization schemes. 

PACS numbers: 11.15Ha, ll.lO.Gh, 12.38.Bx 
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I. INTRODUCTION 



For some time a goal of lattice QCD has been to determine the masses of the quarks. To 
be precise, one would like to quote a value of m(/i), the renormalized mass in the MS scheme, 
at momentum scale /i. This is the convention most often used in the phenomenology of the 
Standard Model and in attempts to treat the Standard Model as the low-energy limit of a 
more fundamental theory. 

In calculations of the hadron spectrum the bare (lattice) mass is a free parameter, which 
is adjusted to match experiment. The MS mass is related to the lattice mass via perturba- 
tion theory. This relation is obtained by computing the quark's pole mass in dimensional 
regularization and in lattice perturbation theory, and then eliminating the pole mass. For 
the light quarks the perturbative matching is well established at the one-loop level. The 
results for the MS mass (in the quenched approximation) form a consistent picture, at least 
when power-law lattice-spacing effects (from the underlying hadron masses) are taken into 
account 

For the charm and bottom quarks lattice artifacts may seem, at first glance, a greater 
worry, because m^a (the quark mass in lattice units) is not necessarily small. This is, 
however, not so. Lattice artifacts take the form 

aSE = a''''bn{mqa){SHn) (1.1) 

where 6Hn is an operator whose matrix elements are (typically) insensitive or mildly sen- 
sitive to the heavy-quark mass, and h{mqa) is a c-number function. In the static [0,^ and 
nonrelativistic [^,^ effective theories, which develop an expansion in l/rriq, Eq. ( |1.1| ) arises 
by design, and the 6„ are bounded for large rriqa. For Wilson-like actions Eq. (|1.1|) also 
holds 0, and the 6„ are bounded for all irtqa. To obtain this result, it is essential to avoid 
expanding around m^a = or around ntq = oo aX every stage of the analysis. 

The implication of Ref. for lattice perturbation theory is that the relationship between 
renormalized and bare quantities is needed for arbitrary m^a. This paper examines mass 
and wave-function renormalization in the class of actions considered in Ref. and gives 
concrete results at the one-loop level for the Wilson and Sheikholeslami-Wohlert 
actions. Results are available in the literature for niqa = P-pT[] and niq ^ oo 0; as 



expected [§[, the new results presented here smoothly connect the two limits. 

The mass dependence of one-loop lattice perturbation theory has been considered before. 
Results from nonrelativistic theories [0,|l3 provide us with cross checks, because (suitable 



combinations of) their results must agree with ours in the static limit, niq oo. At the 
other extreme, terms of order rriqa, from expanding our general results around rriq = 0, 
should recover the results of Sint and Weisz 

This paper is organized as follows: Section |I| discusses the pole in the lattice quark prop- 
agator for all rriqa, to all orders in perturbation theory. The action of Ref. 0, special cases 
of which are the Wilson and Sheikholeslami-Wohlert actions, is reviewed in Sec. The 
main (all orders) results of Sec. H are expanded to first order in in Sec. |V|, and expres- 
sions for one-loop Feynman diagrams are presented. Some notation necessary for simplifying 
numerical evaluation of the one- loop diagrams is given in Sec. [V|. Section |VT| presents numer- 
ical results for the one-loop contributions to the rest mass, the kinetic mass, and the wave 
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function renormalization factor. (Some of the results have appeared previously |[T5|-[T7|.) As 
usual, the dominant contributions come from tadpole diagrams; the results of Sec. ^ are 
improved by "mean-field theory" |T8| in Sec. |V11| . Some technical details are deferred to the 
Appendices. 

The perturbation theory of this paper will be combined with Monte Carlo calculations 
of the quarkonium and heavy-light spectrum, to determine MS masses rhch and rhh, in 
forthcoming publications | |19|| . 



II. RENORMALIZATION TO ALL ORDERS IN PERTURBATION THEORY 

The objective of this section is to derive relations for mass and wave-function renormal- 
ization, for arbitrary values of rriqa. To do so, we shall assume only properties guaranteed 
in the lattice theory, and we shall not assume that the self energy is small. By assuming 
less, we obtain more: our derivation succeeds not only for arbitrary mass, but to all orders 
in perturbation theory as well. 

An outline of our analysis is as follows: We start by anticipating the physical content 
of the quark propagator, viewed as a function of three-momentum and time. This provides 
a template from which one can read off the energy, as a function of three-momentum, and 
the wave-function renormalization factor. We next write down a description of the free 
propagator and the self energy, as functions of the four-momentum, constrained only by 
symmetry and periodicity. The description applies to all lattice theories with Wilson's 
discretization in time |^. This includes the Wilson action |^ (of course), the improvements 
of Ref. PI and of Ref. 0, and some nonrelativistic actions. Then we Fourier transform 
the full propagator from four-momentum to time and three-momentum. The result is a 
sum of terms, one for each pole in the momentum-space propagator. Finally, we focus on 
the pole corresponding to the one-quark state, and read off the energy and wave-function 
renormalization factor from the template anticipated at the outset. 

In our analysis it is unnecessary to assume that the three-momentum, the mass, or even 
the self energy itself, is small. Indeed, the notion of a perturbation arises only to separate the 
inverse full propagator into a free part plus a self energy and to identify a one-quark state in 
the interacting theory. Thus, our results for the pole are valid not only for all masses but also 
to all orders in the gauge interaction. (They are not nonperturbative, because quark and 
gluon states have meaning, in nonabelian gauge theories, only within perturbation theory.) 
Specializing to small three- momenta (in lattice units), we obtain the three main results of 
this section: all-orders formulae for the rest mass (Sec. |TBD, the kinet ic mass (Sec. |II Q ) , and 
the wave-function renormalization factor (Sec. [11 D| ). In subsequent sections, we specialize 
further, first to expressions for the one- loop self energy for the action given in Ref. , and 
later to numerical results for the Wilson and Sheikholeslami-Wohlert actions. 



A. The Quark Pole 

Because of confinement, the true states of QCD are hadrons, not quarks and gluons. In 
perturbation theory, however, one may pretend that quark and gluon states exist. Although 
the aim of this section is to relate the bare mass of lattice QCD to the perturbative pole 
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mass, one should always view the pole mass as an intermediate step. In a final application, 
the pole mass should be related to another regulator mass, such as the MS mass, or to a 
genuine (hadronic) observable. 

Mindful of the preceding caution, we proceed as if quarks and gluons are physical states. 
The starting point is the quark two-point correlation function 

m,p')mp)) = {2nf5ip-p')Git,p), (2.1) 

which defines G{t,p). The fields are those appearing in the functional integral: they are 
bare fields, in a fixed gauge chosen so that G does not vanish trivially. The field ilj{0,p) can 
create from the vacuum not just the one-quark state, but also states with extra gluons or 
extra qq pairs. One thus anticipates 

G(t,p)=Z2(p)e-^(P)l*lQ(p)+/ ^Z,,(p,fc)e-^-(P''=)l*l + ---, (2.2) 

where E{p) denotes the energy of the one-quark state with momentum p, and Egg the energy 
of states with a quark and a gluon. (The quark-gluon states and the multi-particle states 
denoted by the ellipsis will not concern us further.) If the 7 matrix Q{p) is normalized 
according to the condition given in Sec. |IID| , then the residue Z2{p) is the square of the 
amplitude for the bare field to create a physical one-quark state. 

With a Euclidean invariant cutoff the energy would satisfy E{p) = y/p'^ + m^, where m 
is the quark's "pole" mass. With a lattice cutoff, on the other hand, the mass shell is 
distorted. To describe the distorted pole position in a systematic way, one can define a rest 
mass 



Mi=E(0), (2.3) 

a kinetic mass 

M,= f?|V', (2.4) 



dpi 



p=0 



and so on. In general Mi 7^ M2, though as a ^ one should find M2 —>■ Mi. Alternatively, 
at nonzero lattice spacing one can impose M2 = Mi as a requirement on an improved 
action 0. 

Similarly, with a Euclidean invariant cutoff the residue Z2{p) would be a function of 
only, evaluated on shell at = —m?. It is thus a constant independent of p; it is the 
wave-function renormalization factor Z2. With a lattice cutoff, however, the p dependence 
does not drop out of the residue, even on shell. A reasonable definition of the wave- function 
renormalization factor is the residue at vanishing three-momentum 

Z2 = ^2(0). (2.5) 

Then Z2 ip creates the one-quark state with conventional (unit) normalization, at least 
for momenta much lower than the ultraviolet cutoff. 
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To obtain perturbative expressions for E{p) and (p), one starts in momentum space. 
The inverse full propagator is written 



G-\p) = G,'ip)-j:ip), (2.6) 

where the self energy S(p) is the sum of all one-particle irreducible graphs. Given Gq and E 
one obtains the Fourier transform 

Git,p)= r''' ^e^PotG(^po,p) (2.7) 

and compares with Eq. ( p.2|) to obtain E{p), Z2{p), and Q{p)- 

We now introduce suitably general expressions for the propagator. For the lattice theories 
under consideration one can write the inverse free propagator as 

aG^\p)=tI/C{p) + L{p), (2.8) 

where 

KoiPo,p) = sin(poa), 

K{po,p) = K{p), (2.9) 
L{Po,P) = Kp)- cos(poa)- 

The exhibited po dependence corresponds to Wilson's discretization, but the functions K{p) 
and /i(p) depend on the action.[] We assume the action conserves parity, hence K{p) is an 
odd function of p, and ii{p) even. The couplings are not explicit here, but reappear below as 
coefficients in the Taylor expansion of K and around p = 0. For example, /i(0) = 1 + moa, 
where iriQa is the bare mass. 

We decompose the self energy into 7 matrices similarly to Eq. ( |2.8| ) 

aL{p) = I J2 7p^p(p) sin(p,a) + C{p). (2.10) 

p 

With a Euclidean invariant cutoff Ap = A would be a single function for all p, and A 
and G would depend on p"^ only. With a lattice cutoff, however, they are constrained only 
by (hyper)cubic symmetry. For example, symmetry under parity implies that Ap{pQ,p) 
and G{po,p) are even functions of po and p; symmetry under cubic rotations implies that 
Ai(0,p, 0,0) = ^2(0,0,^,0); etc. Furthermore, Ap and G are periodic functions of po, 
with period 27r/a. Below we make no assumptions about the self energy, except for these 
symmetry and periodicity properties. 

Substituting Eqs. ( [^.8| ) and ( [^.10| ) into Eq. ( [^.7| ), and adopting lattice units (a = 1), one 
finds 

G(t,p) = r ^ e^'^'Nip^p) 

' J-n 2-71 2^ {po,p) [cosh £{pq,p) — cospo]' 



For nonrelativistic actions, one takes K = and projects out only the ^(1 + 7o) component. 
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where i^ipo,p) = jJ'ip) - C{po,p), 

N{po,p) = z7osmpo[l - ^o(Po,p)] +iJ2'yi[^i ~ smpiAi{po,p)] - [z/(po,p) - cospo], 

i 

(2.12) 

and 

2i^cosh£:(po,p) = 1 + z/^ + YllKi - Ai{po,p) sinpi]^ 

i 

-sm2po{l-[l-^(po,p)]'}. (2.13) 

For t 7^ one can integrate over pq by changing variables to 2; = e*^"'^^'^*; then one has 
contour integration around the unit circle, and the integral is obtained through the residue 
theorem. 

The integrand has a pole at z = e~^, whenever E solves the implicit equation 

coshE = cosh^(iE,p). (2.14) 

No more compact, general expression for E exists but to set po = iE in Eq. (|2.13|) . (At 
fixed order in perturbation theory one solves Eqs. ( p.l4| ) and ( |2.13| ) iteratively.) Solutions 
of Eq. (|2.14|) are parametrized by the three- momentum p and will be denoted E{p). 

In general, the integrand of Eq. ( |2.11| ) has several poles. In perturbation theory one 
assumes, however, that the self energy is a "small" correction. Then the pole corresponding 
to the quark state must have a residue that does not vanish as Ap, C — * 0. Poles corre- 
sponding to multi-particle states, on the other hand, must have residues that do vanish in 
the absence of an interaction. 

Given an energy satisfying Eq. ( p.l4| ), one expands cosh. S{—i In z,p) in z around 

cosh.S{—i\nz,p) = cosh E + z~^{z — e~^)S{iE,p) sinhE (2.15) 
+ 0{{z-e-y), 

where 

£ = z-^£{-i\nz,p) = -^{po,p). (2.16) 

az I apo 

In applying the residue theorem, the quadratic and higher-order terms drop out, and 
Eq. ( |2.11|) becomes 

^ ^ [1 - Ao{iE,p)][-fo sgntsinh E - ij ■ P{p)] + iy{iE,p) -coshE ^^^^t^ 

2u{iE,p)[l + S{iE,p)]smhE ^ 

+ other residues, 
where, for brevity, E = E{p) and 

, , Kiip) — Aj(iE,p) sin Pi 
1 - Ao{iE,p) 

The chosen pole has a residue that remains when the interaction is turned off; thus, it 
corresponds to the one-quark state. The "other residues" correspond to states other than 
the one-quark state and are disregarded from now on. 
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B. Rest Mass Mi 

To obtain an expression for the rest mass, one sets po = T'E{p) in Eq. (|2.13|) and solves 
for V. One finds 

u{iE, p) = flip) - C{iE, p) 

= coshE + (1 - Ao)\/smh^ E - P^. (2.19) 

Setting p = in Eq. (|2.19| ) yields the implicit equation 

e^^i = l + mo + AoiiMi, 0) sinhMi - C{iMi, 0), (2.20) 

where the parameter 

mo = /i(0) - 1 (2.21) 

is the bare mass. Equation (^1201) , expressing the rest mass to all orders in perturbation 
theory, is the first main result of this section. 

For a massless fermion, the rest mass Mi should vanish. The critical bare mass, which 
induces Mi = 0, is 

moc = C{0,O;moc). (2.22) 

The self energy depends on mo as a parameter, denoted here by the third argument of C. 
Since the lattice actions under consideration do not maintain explicit chiral symmetry, one 
expects moc 7^ 0. When applying Eq. ( |2.20| ) it is useful to take care of this term once and 
for all and write 

e^i = 1 + Mo + Ao{iMi, 0) sinh Mi - C{iMi, 0) (2.23) 
where C{iMi, 0; mo) = C{iMi, 0; mo) — moc, and 

Mo = mo - moc = 7^ - (2.24) 

In practice, one determines moc (or Kc) nonperturbatively in Monte Carlo calculations and 
treats Mq (rather than mo) independently of Qq in perturbation theory. 



C. Kinetic Mass M2 

From its definition [Eq. (|2.4| )] the kinetic mass requires two derivatives with respect to pi. 
Because the derivatives are applied to on-shell self-energy functions, the total derivative with 
respect to pi includes an explicit part and an implicit part through the dependence on E, 

d d dE d , , 

l- = 7^ + 'f^l^- 2.25 
api opi opi opo 



Differentiating Eq. ( p. 191) twice yields 
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e^^ - Ao(zMi, 0) coshMi _ ^ ^ , , [C - A,{iM,, 0)f 

where 



= ^ [Ao{iEip),p) sinhMi - CiiEip),p)] . (2.27) 
"Pi 



The quantities 



C^K[iO), (2.28) 
r,C = /i"(0), (2.29) 

with primes denoting differentiation with respect to pi, are couphngs parametrizing the 
action. Equation ( pj.26| ), expressing the kinetic mass to all orders in perturbation theory, is 
the second main result of this section. 

Without the interaction that generates the self energy, one can solve Eq. (|2.26| ) in closed 
form: M^"' = m2(M|'''), where 

, , e™ sinh m 
m2im) = . 2.30 

C + TgQ smhm 



This expression suggests defining a kinetic-mass renormalization factor Z^^ through 

M2 



'M2 



m2(Mi 



(2.31) 



which captures the radiative corrections to M2 not shared by Mi.0 In perturbation theory a 
quark's kinetic and rest masses are on-shell observables. By definition the ratio is also 
an on-shell quantity. It is, therefore, a useful diagnostic of cutoff effects in the continuum 
limit. Indeed, 

lim Zmo = 1 (2.32) 

to all orders in perturbation theory. 

For Mia 7^ the rest and kinetic masses are not, as a rule, equal. To construct a mass- 
dependent improved action, as in Ref. 0, one sets M2 = Mi in Eq. ( p.26|) and solves for 



a condition on the coupling parametrized by r^. In the following we do not, however, 
assume that such conditions have been imposed. 



^Please note that it is the ah-orders rest mass that appears as the argument of the function m2- 
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D. Wave-function Renormalization 



Comparing Eq. ( p.lTD with Eq. (p.2|) one identifies the 7 matrix function of p that 
multiphes e"'^*-^'"*' as Z2Q. To fix the normahzation of Q, first write Q = sgnt'joQQ — 
i-y ■ Q + R. (Note that Qg ~ — i?^ = on shell.) In the canonical normalization Qq = \. 
Thus, 

= . (2.33) 

v{iE,p)[l + £{iE,p)] 



A more explicit expression may be obtained by using Eqs. ( p.l3| ) and ( |2.14| ) to eliminate 
fiiE^p) and p£{iE,p). One finds 



Z2{p)-^ = (1 -Ao)cosh^+ \/smh^ E - + AosinhE (2.34) 

AjPj sin pj 
sinh E ' 



C^Jl- Pysinh^ E -Y, 



with all self-energy functions evaluated at {iE{p),p). Setting p = to obtain the wave- 
function renormalization factor, as discussed above, one finds 

= z^ioy' 

= e*^i - Aq cosh Ml + io sinh Mi - C, (2.35) 

where all self-energy functions are evaluated at (iMi,0). Note that every term on the 
right-hand side of Eq. ( p.35 ) is of order e*^^ in the large mass limit, just as at tree level. 



Equation (|2.35|), expressing the wave-function renormalization factor for arbitrary values 



of Mia, to all orders in perturbation theory, is the third, and final, main result of this 
section. 

In most gauges, the wave-function renormalization factor is infrared divergent. The 
divergence cancels against vertex renormalization factors, when a physical combination, 
such as the full renormalization of a current, is considered. 



III. THE LATTICE ACTION 

We consider the action S = Sq + Sb + Se of Ref. 0, namely 

^0 = moa' H^M^) + y E V^(^) [(1 + 7o)^o" - (1 - lo)D+]ij{x) 

X X 

+ Ca^ E '^(^h ■ Dipix) - IvsCa^ ^ Hx)A''^^ip{x), (3.1) 

X X 

where the covariant difference operators are 

aD^ip{x) = ±[U±o{x)'ilj{x ± aO) - ^(x)], (3.2) 



aDiip{x) = ^[Uiip{x + ai) — U^iip^x — a?)], (3.3) 
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and 



a^A^^^ipix) = J2 iPii^i^ + + U^i-^ix - at) - 2iIj{x)]. (3.4) 

i 

The action 5*0 has cutoff artifacts of order a, which can be canceled by the interactions 

SB = -^a''cBCj2^i^)^-Bix)^ix), (3.5) 

X 

Se = -la^csC J2 ■ E{x)^Ij{x), (3.6) 

X 

for appropriate adjustments of and ce- The chromomagnetic and chromoelectric fields 
are given in Ref. [^]. 

Special cases of this action are the Wilson action which sets Vs = ( = cb = ce = 0; 
and the Sheikholeslami-Wohlert action which sets = ( = 1, cb = ce = csw- But to 
remove lattice artifacts for arbitrary masses, the couplings r^, (, cb and ce must be taken 
to depend on nioa For the purposes of this paper, however, the additional couplings are 
taken as free parameters. 

We note here the elements of the free propagator introduced in Eq. ( p.9| ). From Eq. ( p.l|) 
one finds 

Ki{p) = (smpitt, (3.7) 
fi{p) = 1 + rrioa + \rsCp^a?, (3.8) 

where Pi = (2/a) sin |pja. Thus, the notation for the couplings coincides with that in 
Eqs. ([23l| ), (PI), and 



As the mass tends to infinity, all actions described by Sq + Sb + Se lead, up to an 
unphysical factor, to the same quark propagator — a Wilson line. Perturbative corrections 
to masses and vertices must respect this universal static limit, and, therefore, they must 
tend to a universal value. This limiting behavior is a helpful check. 

IV. THE SELF ENERGY TO ONE LOOP 

The analysis of Sec. is valid to all orders in perturbation theory. We now develop 
expansions in for the main results, concentrating on the one-loop approximation. We also 
present our expressions for the one-loop self energy. 

A. Perturbative Series 

In perturbation theory the self energy is expanded 

oo 

S(p)=E^?o'sW(p), (4.1) 
1=1 
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and similarly for the functions Ap and C. As a consequence, the rest mass has an expansion 



V] 



1=0 



where the tree level m|°^ = log(l + Mq). The one-loop coefficient, from Eq. (|2.2CI| ), is 



In the massless limit 



A™ (iMf , 0) sinh Mf - C^^^ (iMf , 0) 



-Al- 



io] 



[1] 



M 



.[1] _ 

^0 



dniQ 



(4.2) 



(4.3) 



(4.4) 



mo=0 



which is the same result as in massless derivations 0. In the static limit, Aq tends to a 
constant and C to e*^^ times a constant. Thus, m|^'(cx)) is finite. Moreover, it is the same 
for all actions under consideration. 

The kinetic-mass renormalization factor has an expansion 



'M2 



21 

M2 



(4.5) 



1=1 



The one-loop coefficient, from Eq. ( p.26|) , is 

2CAS'l(iMf , 0) - C^A\^\iMf'\0) - DW(0) sinhMf 



.[1] 



C + rX sinh Mj 



In the massless limit, vanishes at least as fast as (Ma)^ In(Ma). In the static limit 



[0] 



-A\l\iMf\0) cosh Mr e 



[0] 



(4.6) 



'Ah 



'M2 

—D^^^/{rsC) — \a\^\ again finite. Moreover, it too is the same for all actions under 
consideration, and it can be compared to the same combination in nonrelativistic QCD. 

At tree level many have noticed that the quark propagator's residue ^2(0) = e~*^^. This 
dominant (large) mass dependence persists in individual loop diagrams and, as shown in 
Eq. ( p.35| ), to all orders. To isolate the subleading mass dependence of the wave-function 
renormalization, we develop the expansion as follows: 



(4.7) 



1=1 



with the physical (all-orders) rest mass in the exponent on the left-hand side. We apologize 
for a notation in which is not the perturbative coefficient of Z2, but this way the 
have only mild mass dependence. The one-loop coefficient, from Eq. ( p.35|) , is 



'[1] 



(iMf , 0) cosh Mf - ij)'^ (iMf , 0) sinh Mf ' + C^'^ (iMf ' , 0) 



-M- 



[0] 



(4.8) 



In the massless limit, the right-hand side reduces to the well-known result Aq — mAq + (7^^' . 
(As M ^ 0, const.) In the static limit, each term on the right-hand side of 

approaches a (universal) value, and thus the sum Z^^' does too. 
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(a) 



p+k 



(b) 



FIG. 1. Feynman diagrams for the one-loop self-energy: (a) rainbow, (b) tadpole. 



B. One-Loop Diagrams 



It is straightforward to derive Feynman rules for Sq + Sb + Se- They are listed in 
Appendix ^ To one loop the self energy is given by the Feynman diagrams in Fig. |l[ We 
use the Wilson gauge action [Q. The rainbow diagram, Fig. |I]a, represents the contribution 



F{p, k) 



i2TT)^ k^[K^{p + k) + L^{p + k)y 



(4.9) 



with Kp and L specified by Eqs. (|2.9|), (^.7]), and (P78|). The numerator F{p, k), which follows 
from 7-matrix algebra, is given in Appendix |B[ The color factor Cp = {N^ — l)/{2Nc) 
[= 4/3 for SU(3)]. The tadpole diagram. Fig. P^, represents a much simpler contribution; 
in Feynman gauge 



which is independent of and ce- 

On shell the total one-loop self energy is gauge independent 



d'^k 1 
(2^F 



(4.10) 



Therefore, the one-loop 

radiative corrections mJ;^' and Z^^^ are gauge independent, as one expects. Derivatives of 
the (off-shell) self energy do, however, depend on the gauge parameter. Therefore, the 
wave-function renormalization factor does depend on the gauge; below we present the re- 
sult in Feynman gauge. Note, furthermore, that terms arising from Sb and Se are gauge 
independent (at the one-loop level). 

By infrared power counting, one expects the one- loop self energy S^^^ to contain logarith- 
mic nonanalyticity as Ma — > 0. The leading nonanalyticity is the same as in a Pauli-Villars 
regulator. The latter amounts to using the gluon propagator 



1 



(4.11) 



where the gluon mass A serves as an infrared regulator. As with the lattice calculation, we 
do not necessarily assume that Ma is small, and we usually set a = 1. Thus, 
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As in Eq. ( p.lO[ ) we write 



^PV 



167r2 {p + ky 



ipv -t- '-^pv 



Below it is sometimes convenient to set 



Cpv(P^; m) = mB^pllp"^] m), 



(4.12) 



(4.13) 



(4.14) 



because the Pauli-Villars regulator does not break chiral symmetry. Equation ( [4.1 Ij ) specifies 
a Euclidean-invariant cutoff, so Ap^ and Bpy are functions of p^, rather than of {po,p)- 
Logarithms arise in Eq. (gl2D from the region /c^ <^ 1. 

Below we exploit the similarities of the two regulators to isolate analytically terms of the 
form In Ma and Ma In Ma. 



V. SUMMARY OF NUMERICAL METHODS 

In this section we outline the numerical procedures used to evaluate the loop integrals. 
In particular, we isolate from Mi'\ Z^lj^, and ^2^^ parts that are easy to compute numerically. 
The notation introduced is needed below to obtain the one-loop coefficients, as a function 
of mass, from the tables in Sec. [V^. Some other technical details are deferred to appendices. 

The remainder of this paper focuses on one-loop renormalization for the action used in 
Monte Carlo calculations, namely Vg = C = 1 and cb = ce = csw- For future flexibility it is 
useful to classify the results as a (second-order) polynomial in csw, for example 

AW=AW^ncswAW«+4w4^1(^^ (5.1) 



A. Numerical Integration 

Because one must anal5d;ically continue the self energy from real po to imaginary iE{p), 
it proved wise to carry out the integration over ko analytically. A full discussion of this 
technicality is in Appendix 0. Here we focus on the remaining integration over k. 

In numerical evaluation of the lattice integrals, it is helpful to compute the difference 
between the lattice and Pauli-Villars regulated integrals. The key is to subtract the two 
integrands in momentum space and to add the analytical expression for the Pauli-Villars 
integral afterwards. Then the numerical integration package does not need to uncover the 
logarithmic singularities. At small k the lattice integrands take the form 

"^'""^ (5.2) 



2sinhMr|A;| 



The denominator of the Pauli-Villars integrand takes the same form if one sets m = sinh m|°' 
in Eq. ( |4.12| ). If one also multiplies the Pauli-Villars integrand by the function p, the 
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numerical integration package has an even easier job, because then the subtraction removes 
contributions of the form Main Ma (for small Ma). 
Thus, let 



"0 



4[l](n) 
^0 



4l](n)^^[l](n)_p(n)(^[0])^[l]^ 



(5.3) 
(5.4) 
(5.5) 



The Pauli-Villars subtractions on the right-hand sides are done on the integrands. They are 
needed for the Wilson-action (n = 0) contribution only. We find 



pS(M) = |(3e-^^-e^^) 



pJ)(M) = e-*'(l - sinhM) 



pP{M) = U3e-'^^ + e 



(5.6) 
(5.7) 
(5.8) 



and pj^^ = otherwise. The subtraction C{iMi;mQ) — C(0;moc) = C* is also done on the 
integrand. 

The kinetic mass requires also the function D, defined in Eq. (|2.27|) . The total derivative 
with respect to pi acts on the energy E{p) and the explicit pi-dependence. Each generates 
a severely divergent peak in the integrand of Z}[^](°) at somewhat different infrared locations. 
After integrating, the infrared divergences cancel exactly. To make numerical integration 
easier, it is again prudent to let 



(5.9) 



to cancel the peaks against integrands from the Pauli-Villars regulator. The function p£, 
is not needed below, however, because the integral Dpy vanishes identically, but we take 
Pd{M) = e^7[coshM + |(e-^^ - 1) sinhM]^. 



For the wave-function renormalization factor we compute 



Win) 



coshMj'^ 



Ji]W 



sinh MP + d 



WW 



[0] 



(5.10) 



where do^"-"'^ and cl^l^") are defined analogously to ag^'*-"^ and c^^^ The Z2'^^"'^ are finite 
as A — > and as Ma —^0. In any physical quantity the wave-function renormalization is 
combined with vertex corrections in an infrared-finite way. Thus, 2:2^^ is a suitable synopsis 
of the lattice renormalization. 

The subtractions permit a numerical evaluation of Og^' a^^^ ct^l*^"\ d^^^^"'\ and ^2^' 
with gluon mass A = 0. With the subtracted integrals in hand, the lattice self energy can 
be reconstructed with the closed forms for Apy, Bpy, Apy, and Bpy, given in Appendix 0. 



B. Chebyshev Approximation 

With the adaptive integration routine VEGAS we evaluate aQ^^"'\ a^i^^"\ c^'^^'^"'\ d^^^''"'\ 
and Z2^^"^ at 51 values of the mass, chosen such that 
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tanhMf' = |(1 + x^), k = 0, 50, 



(5.11) 



where the Xk = cos[7r(/c + |)/51] are the zeroes of the 51st Chebyshev polynomial. This 
procedure allows us to combine the individual evaluations into a Chebyshev approximation 
to the exact result. Let 



N-l 



(5.12) 



fc=0 



where = 51 and Tj{x) = cos(jcos x) is the jth Chebyshev polynomial. Then 



m—l 



fix) ~ i/o + E fkTkix 



k=l 



(5.13) 



is (expected to be) a good approximation even for m < N . Since |Tj(x)| < 1 the utility of the 
approximation can be ascertained from inspecting the fj. Section ^ gives tables with the 
first several coefficients of Chebyshev expansions. All 51 values of f{xk) and fj are available 
on the Worldwide Web at [http : / /www-theory . f nal . gov/people/ ask/ self -energy/ . 
The lattice self-energy functions have also been obtained over a wide range of masses |T7 



without the Pauli-Villars subtractions. The results agree, of course, but with the subtrac- 
tions one can reach better precision more quickly. 

VI. ONE-LOOP RESULTS 
A. Critical Bare Mass moc 



For completeness we give here our result for the one-loop bare mass that makes the 
physical masses vanish: 



[1] 

Oc 



CW(0, 0) = Cf [-0.325714(5) + 0.086964(9)csw + 0.036190(2)^^ 



(6.1) 



Errors on the least significant digit (s) from numerical integration are given in parentheses. 



Equation (|6.1| ) agrees with published values: —C^O. 325789(3) at csw = 00 and — Ci?0.20(l) 
at Csw = 1 Hoi • The individual coefficients of csw agree with Ref . |]lT | . 



B. Rest Mass Mi 

Figure § shows the one- loop correction to the rest mass Mi. We present results for three 
values of csw^ (Wilson action), 1 (tree-level improvement), and 1.4 (a typical mean-field 
estimate of csw)- As expected, M^^ smoothly connects to the massless and static limits. 
As Ml oo all curves approach the same limiting value, Mf^(oo) = Cf0.1261(2), which 
agrees with the value Ci^O. 1263(1) obtained directly in the static limit We are able to 
reproduce the result of Ref. [^, which considers only the Wilson action, if we omit the 
tadpole diagram's contribution e[|^1. 
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0.2 



I I I I I 




1 2 3 4 5 

Ml°] = ln(l +Mo) 



FIG. 2. Plot of M^^ vs. Mf' for 

csw = (circles), 1 (triangles), and 1.4 (squares). The 
diamond marks the static limit To illustrate that the overall shape is not an artifact of lattice 
field theory, the curve shows M]L^. Here Cp = 4/3. 

For Ml a < 1 the additive form is not illuminating, because there m|^^ oc Mf\ It is 
convenient to define the rest-mass renormalization factoif] 

Mia 



tanhMf a 

To build the one-loop renormalization factor from the VEGAS integrals let 

[l](n) ^ Mfl / [l](n) 



and then 



e-<' (aW(") cosh Mfl - cl^^^") cothM^ 



'Ml 



"Ml 



+ e-^^i~' coshMr p^2l{Mf^)^VY - P?(Mf )4'i 



(6.2) 



(6.3) 



(6.4) 



Figure ^ shows the z^l}^^'^^ and Table | contains the first 15 coefficients of their Chebyshev 
expansions. 

In the massless limit there are several checks in the literature. As Mia — > 0, we find 



0.10726(15) + 0.04901(2)csw - 0.008735(5)4^ - ttt^ ln(Mf a 



(6.5) 



Note the appearance of the logarithm, multiplied by the one-loop anomalous dimension. The 
finite part agrees well with pubhshed values: (7^0.107347(5) at csw = 0, and Cj^O. 1474(4) 
at Csw = 1 IM- The individual coefficients of csw agree with Ref. [111 . 



^The denominator tanhMj is handy because tanhm = m + 0{m^) at small m, yet tanhoo = 1. 
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I I I I I I I I I I I I I I I I I I I I I I I 
1 2 3 4 5 

M^l = ln(l + Mo) 

FIG. 3. Plot of VS. Mj°^ as coefficients of csw- The curves indicate the fifteen-term Cheby- 
shev interpolation. Note that ^j^^^^"* (the coefficient of Cg^) is multiplied by ten. Here Cp = 4/3. 



TABLE I. Coefficients of Chebyshev polynomials Tj{x) for z^^l''' . Here Cp = 4/3. 



j 


^[11(0) 









0.222802 


0.0868169 


-0.00278643 


1 


-0.0117152 


-0.0300771 


0.00744513 


2 


0.00883238 


-0.0052349 


-0.00297984 


3 


0.00216047 


-0.000493063 


-0.00039969 


4 


0.00141579 


-0.00245913 


-0.000478208 


5 


0.000994781 


-0.000371327 


-0.000358493 


6 


0.000590983 


-0.000965165 


-0.000237546 


7 


0.000467076 


-0.000311904 


-0.000198283 


8 


0.000321860 


-0.000497263 


-0.000147608 


9 


0.000256113 


-0.000225545 


-0.000123028 


10 


0.000195137 


-0.000302366 


-9.86053x10-5 


11 


0.000159702 


-0.000161754 


-8.34899x10-5 


12 


0.000125988 


-0.000203734 


-6.96978x10-5 


13 


0.000106567 


-0.000120868 


-5.99954x10-5 


14 


8.55174x10-5 


-0.000146607 


-5.15038x10-5 


15 


7.61065x10-5 


-9.24268x10-5 


-4.50976x10-5 
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Recently Sint and Weisz |14] have computed the next term in the expansion of Z^^j (Mia) 
around Mia = 0, for csw = 1- They find the coefficient of Mia to be —CirO. 07217(2). 
Fitting our resuhs for Mia < 0.1, we find —(7^0.0720(7), less precise, but in agree- 
ment. Reference fl^ does not report a contribution of order MiolnMia; for csw = 1 



one would expect it to drop out. Our Pauli-Villars subtractions isolate from z]^}^^^ a contri- 
bution CF[6/167r^]MialnMia, and our fits find nothing more of order Mialn Mia in Z^^}^^\ 
But Z^^jj^^ contains precisely the same contribution with the opposite sign, so the total drops 



out when csw = 1- This exercise verifies, as in Refs. [1C,11], that the tree-level improved 
action removes terms of order a In a. 

These checks are reassuring, but the main result is the full mass dependence, embodied 
in Figs. and ^ and in Table |. To proceed from our numerical results to z]lj^: 

1. reconstitute adequate approximations to the z^l}^^'' from Table |I|; 



2. evaluate them at the desired value of tanhMj ; 

3. accumulate the polynomial in csw; 

4. add e"^''^' ' cosh m|°' (paq^pv ~ PcBpy) to restore the Pauli-Villars subtraction. 
The full one-loop approximation to the rest mass is then 

Ml = m1°1 + ^^4/1 tanh mJ°^ . (6.6) 

In a straightforward application of bare perturbation theory the expansion parameter g"^ 
would be the bare coupling g^. It is possible, however, to choose a better expansion param- 
eter Further discussion of this issue will appear elsewhere. 



C. Kinetic Mass M2 

Figure § shows the one- loop renormalization of the kinetic mass for csw = and 1. (The 
variation with csw is too weak to distinguish 1.4 from 1.) Again, Z^^^ smoothly connects 
to the massless and static limits. The separate coefficients of csw are plotted in Fig. |, and 
their first fifteen Chebyshev coefficients are listed in Table |T|. 

In the static limit we find ^/^(oo) = -Cj.0.0745(l) = -0.0993(1), which agrees with 
—0.0998(4), the value in nonrelativistic QCD |ll2| , p^ for our definition of Zm2- In the massless 
limit we verify Z^lj^{0) = 0, because there M2 = Mi. 

These checks are again reassuring, but the main result is the full mass dependence, 
embodied in Figs. ^ and ^ and in Table |I|. To proceed from our numerical results to z]lj^: 

1. reconstitute adequate approximations to the Z^^l^^ from Table 

2. evaluate them at the desired tanh m|°^; 

3. accumulate the polynomial in csw- 
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-0.1 ^ 



2 3 4 

mP = Iog(l+Mo) 



FIG. 4. Plot of Zjj^ vs. a for csw 
static limit pi. Here Cp = 4/3. 



(circles) and 1 (triangles). The diamond marks the 




12 3 4 

mI,°J = log(l+Mo) 



FIG. 5. Plot of ZJ5^^"^ vs. a. The curves indicate the fifteen-term Chebyshev interpolation. 
Note that the coefficients of and Cg^v ^^e multiplied by ten. Here Cp = 4/3. 
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TABLE II. Coefficients of Chebyshev 


polynomials TAx) for zW^"'^ 


Here Cp = 4/3. 


3 











-0.0631437 


-0.00518218 


-0.000697017 


1 


-0.0409363 


-0.00266278 


-0.000293363 


2 


-0.0113513 


0.000649031 


0.000188227 


3 


-0.00411835 


0.00104410 


0.000188092 


4 


-0.00293037 


0.000578677 


9.36398x10-5 


5 


-0.00157589 


0.000435851 


5.85129x10-5 


6 


-0.00118956 


0.000331139 


3.40842x10-5 


7 


-0.00081670 


0.000251324 


2.23707x10-5 


8 


-0.000630783 


0.000204140 


1.47425x10-5 


9 


-0.000490657 


0.000162222 


1.02279x10-5 


10 


-0.000396438 


0.000137218 


0.71930x10-5 


11 


-0.000319830 


0.000112310 


0.52123x10-5 


12 


-0.000273443 


9.74213x10-5 


0.39003x10-5 


13 


-0.000223656 


8.30367x10-5 


0.29204x10-5 


14 


-0.000198908 


7.12494x10-5 


0.21814x10-5 


15 


-0.000163911 


6.31163x10-5 


0.17610x10-5 



The full one-loop approximation to the kinetic mass is then 

M2 = m2(Mf + glM?^) (l + g^Z^^^) . (6.7) 

Again, it may be appropriate to choose optimal expansion parameters gf and as will be 
discussed elsewhere. 



D. Wave-function Renormalization Z2 

Because the wave-function renormalization factor's full correction ^2^' has an infrared 
divergence for all values of the quark mass, we present results for the subtracted form, as 
defined in Eq. ( |5.1CI| ). Figure ^ shows the one- loop correction Z2'\ in Feynman gauge, for 
csw = 0, 1, and 1.4. Once again, Zg^' smoothly connects the massless and static limits. The 
separate coefficients of csw are plotted in Fig. |^, and their first fifteen Chebyshev coefficients 
are listed in Table fT|. 

In the static limit we find (00) + C^p/ievr^] In = C^O. 1548(5) = 0.2064(7), which 
agrees with 0.2067(1), the value of Refs. [0,|12| for our definition of Z2. 
In the massless limit we find 

= C 



0.05608(7) - 0.014239(6)csw - 0.008844(3)4^ (6-8) 
+ (31n(Mf a)2-21n(Aa)2) /167r2" . 

Note the appearance of a logarithm of as well as the infrared divergence. The finite 
part agrees well with published values: 6*^0.056057(2) at csw = H] and 6*^0.0329(3) 
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0.25 




Q I I 

1 2 3 4 5 

mIi°J = log(l+Mo) 



FIG. 6. Plot of ^2 vs. M[ for csw = (circles), 1 (triangles), and 1.4 (squares). The diamond 
marks the static limit ||3|. Here Cp = 4/3. 




1 2 3 4 5 

mI'^J = log(l+Mo) 



FIG. 7. Plot of 4 ^^^ VS. Mj'^Ja. The curves indicate the fifteen-term Chebyshev interpolation. 
Note that the coefficients of c^^^ and Cg^ are multiplied by ten. Here Cp = 4/3. 
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TABLE III. Coefficients of Chebysliev 


polynomials Tj{x) for 22^^^*^^ 


Here Cp = 4/3. 


3 


ri 1 (r\\ 


^2 


^2 





0.305389 


-0.033722 


-0.0111230 


1 


0.0392549 


0.0108371 


0.00554548 


2 


0.00310868 


0.00597098 


-0.000324921 


3 


0.00393177 


-0.00107846 


0.000324902 


4 


0.00138399 


0.00096398 


-1.79973x10-5 


5 


0.00134333 


-0.000166528 


2.44776x10-5 


6 


0.000666205 


0.000233668 


0.38057x10-5 


7 


0.000655838 


-4.30696x10"^ 


0.283444x10-5 


8 


0.000391408 


8.67030x10-5 


0.156463x10-5 


9 


0.000380178 


-1.79973x10-5 


0.062247x10-5 


10 


0.000255127 


4.08834x10-5 


0.038319x10-5 


11 


0.000252910 


-1.04220x10-5 


0.020393x10-5 


12 


0.000181397 


2.24672x10-5 


0.007974x10-5 


13 


0.000172181 


-0.68263x10-5 


0.005127x10-5 


14 


0.000122340 


1.43152x10-5 


0.001010x10-5 


15 


0.000120061 


-0.53572x10-5 


-0.004476x10-5 



at csw = 1 113 • '^^^^ individual coefficients of csw agree with Ref. |TT|. For comparison 
with Figs. ^ and 0, note that though 2:2^"^*^'* does not contain the logarithms, it is larger 
by Cj.[9/327r2] than the finite part of 

Once again, the checks are reassuring, but the main result is the full mass dependence. 
In practice, the wave-function renormalization factor is used only combined with vertex 
renormalization factors, in ways such that the infrared divergences cancel. When calculating 
a vertex renormalization factor, one should isolate all infrared divergences analytically, as 
we have done here, and then assemble the pieces so that the cancellation is explicit. If 
one chooses an on-shell renormalization scheme, dependence on the gauge parameter should 
cancel as well. 



VII. IMPROVED PERTURBATION THEORY 

The previous section presented results for the (mass-dependent) one-loop coefficients 
in Eqs. (^^), ( |4.5| ), and ( |4.7| ). Truncating the series at one loop, without further ado, is 
unlikely to be a good approximation to the full series, however. The one-loop self energy 
is large, owing to the contribution of the tadpole diagram in Feynman gauge, s'^^-j. This 
feature will persist at higher orders. Similarly, tadpole diagrams make the bare coupling 
much smaller than other, more physical measures of the gauge interaction. Thus, one must 
be wary of these series; they are characterized by an unusually small expansion parameter, 
counteracted by unusually large coefficients. 

These observations suggest rearranging perturbation series so that tadpole diagrams 



cancel each other in final results [18|. The rearrangement is achieved by defining new 
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couplings and (re)normalization factors. Wherever the gauge field appears in the action 
or in operators, substitute 

U,-^Uo[Ujuol (7.1) 

where the new parameter mq, the mean link, is a gauge-invariant average of the link matri- 
ces U^. The perturbative expansion of Uq, like any involving f/^, is dominated by its tadpole 
diagram. In the perturbative expansion of U^j^/uq, on the other hand, the tadpole diagrams 
cancel. The price for arranging the cancellation is that the first factor of Uq in Eq. ( [7.1| ) 
must be evaluated nonperturbatively. 

In the following discussion we shall denote the tadpole-improved gauge coupling as g^, 
without defining it explicitly. The definition of affects higher-order corrections and nu- 
merical evaluation of a series truncated at one loop (or beyond). But it does not affect the 
one-loop self energy, which is the main focus here. On the other hand, tadpole improvement 
of the bare mass and of the field normalization do alter the one-loop coefficients, and they 
are considered below. 



A. Rest Mass Mi 

The rest mass is primarily sensitive to the bare mass, particularly the subtracted version, 
Mq = mo — ruQc- To derive a tadpole- improved bare mass, it is easiest to apply Eq. ( |7. Ij ) to 
the hopping-parameter form of the action. Our action 5*0 has two hopping parameters, 

" 2[l + (rf-l)r,C + mo] ^'^■^^ 

for temporal hops and Kg = (nt for spatial hops. After rescaling the fields by \ph^ti the link 
matrices appear only in the combinations /ttf/o and nJJi. In view of Eq. ( |7. 1|) one defines the 
tadpole- improved hopping parameters kt^^ = WoKt,^ to absorb the first factor of Uq. Undoing 
the rescaling, now with \/2kt, and reassembling the mass form leads to 

Mo = mo- moc = Mq/uq. (7.3) 

The numerical value of the tadpole-improved mass Mq relies on two nonperturbatively de- 
termined parameters, moc and Uq- 

This rearrangement does not alter Eq. ( |2.23| ). Before developing the perturbation series, 
however, one sets Mq = uqMq, treats Mq as the variable independent of g"^, and expands the 
explicit uq as well as the self-energy functions. Thus, 

oo 

M, = Y.f'Mf\ (7.4) 

1=0 



where Mp = ln(l + Mq) and 



1 + Mo 
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1 2 3 4 5 



FIG. 8. Plot comparing z^^j^a vs. Mj^^a (open symbols) and z^^^a vs. Mj^'a (closed symbols). 
The static limit is indicated by diamonds 

Since is negative, the (positive) coefficient m|^' is reduced. Note also that the mass 
dependence of the tadpole improvement connects smoothly to the small- and large-mass 
limits. 

To illustrate the efficacy of tadpole improvement we take the mean link to be uq = 
{UnY^^-i where Un denotes the product of link matrices around a plaquette, = — Cir/16. 
Figure ^ shows explicitly that is significantly smaller than , when uq is defined this 
way. (This choice of Uq does not modify the terms proportional to csw and Cg^.) For csw = 1 
the mass dependence at Ma = is also less steep: -Cir [0.07217 - 1/32] = -CfO.04088. 
Other definitions of uq produce a qualitatively similar reduction. 




B. Kinetic Mass M2 

Because the kinetic-mass renormalization factor is defined as a ratio of on-shell (and 
therefore physical) quantities, it should not be surprising that its one-loop correction Z^l}^ 
is unaffected by the tadpole diagram. More concretely, one can track analytically the con- 
tribution of s'],'^ through Eq. ( |4.6| ), to verify that it drops out (for all Mf^). Thus, there are 

only two changes to improve the one-loop estimate of the kinetic mass: Use Mf^ + g'^Mi^ 
as the argument of m2 in Eq. ( |6.7| ), and use an improved coupling in the expansion of Zm2- 



C. Wave-function Renormalization Z2 

In the rescaling from the hopping-parameter form of the action to the mass form, the 
tadpole- improved version carries an additional factor of Mq, because k = uqk. Thus, the 
renormalization factor of the tadpole-improved wave function is Z2 = U0Z2. Following 
Eq. ( [4.7| ) we develop the perturbative series 
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FIG. 9. Plot comparing ^2 vs. M[ a (open symbols) and ^2 vs. M[ a (closed symbols). 

oo 

e*^^"Z2 = l + X]^'^z{l (7.6) 

1=1 

Comparing the expansions one easily finds 

^f'=4^^+4''. (7.7) 

Figure ^ compares ^2^^ = Zg^'+Wg^' to ^2^^ (in Feynman gauge). Once again, Zg^' is significantly 
smaller. 

VIII. CONCLUSIONS 

The work presented here is the first thorough study of renormalization in the approach of 
Ref. to massive lattice fermions. (Some of our results have appeared earlier ||15|-|17|.) On 
the one hand, the specific, numerical results can be combined with Monte Carlo calculations 
of heavy-quark spectra to determine the heavy quarks' masses. On the other hand, the 
techniques used to arrive at those results can be extended as the renormalization program is 
extended to vertex corrections. The most timely example of the latter is the renormalization 
of (improved) vector and axial- vector currents p2|, which are needed to obtain heavy- light 
meson decay constants and semi-leptonic form factors. 

The theoretical analysis of Sec. || examines renormalization to all orders in the gauge 
coupling and in the fermion mass. Although the focus is on the fermion propagator, the 
analysis serves as a model for the renormalization of currents as well. For currents one 
would introduce appropriate vertex functions, the analogs of the self-energy functions Ap 
and C, again constrained only by periodicity and symmetry properties. Then one would 
Fourier transform the fourth component of each external momentum to obtain the on-shell 
correlation function (for quark states) to all orders in the mass and the gauge coupling. Just 
as here, the all-orders formulae would provide useful insights, and they would be convenient 
for developing expansions to any desired order. 
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In our numerical work we subtract from the self-energy functions corresponding functions 
from a Pauli-Villars regulator. In this way we are able to isolate mass and infrared singu- 
larities, including, for small Ma, contributions of the form Main Ma. The Pauli-Villars 
functions can be written in closed form, so we have analytical control of the singularities. 
Furthermore, by applying the subtraction at the integrand level, numerical integration of 
the nonsingular part is simplified. For example, one can carry out the integration without an 
infrared regulator. These techniques should continue to be helpful for the renormalization 
of currents and other operators. 

The numerical results presented in Sec. |VT| demonstrates that the mass dependence of 
renormalization factors smoothly connects the massless and static limits. This was expected 
from the free theory and general arguments but it is satisfying to make it explicit. Also 
expected, but now explicit, is the result in Sec. |Vll| that tadpole improvement |]18| reduces the 



size of the one-loop coefficients for all masses. In the case of the rest-mass renormalization 
factor, tadpole improvement also makes the mass dependence even smoother, cf. Fig. |^. 

Our results, especially with tadpole improvement, will be useful in determining the MS 
masses of heavy quarks [0. Additionally, such a determination will require at least an 
evaluation of the optimal scale ||18| (in progress), and, to achieve better than 5% accuracy, the 
complete two-loop calculation. Nevertheless, the approach taken here to the renormalization 
of massive lattice gauge theories is a necessary first step. 
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APPENDIX A: FEYNMAN RULES 

This is the first paper to address perturbation theory for the action discussed in Sec. |T|, 
so we give here the propagators and vertices needed for the one-loop self energy. Figure ^ 
defines indices and momentum flow. Then 

Fig. 10a = Go{p)^, = — — (Al) 

«7o smpo + <7 ■ !3[p) + 1 + niQ - cospo + 2^sCp 

Xafe / k k \ 

Fig. [TD|b = AJ(A:) = [6,^ - (1 - a)^j , (A2) 
Fig. [igc = -gotiA'jo cos \{p + p')o - i sin i(p + p)o + hcECcror cos ^ko sin kr], (A3) 
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Fig. |TO|d = -goCtij[jm cos i(p + p')m - ivs sin i(p + p')^ 

+ \cB(Tmn COS \km sill + |cE(7mO COS \kra sill /Cq] 



^oH^"' [cos i(p + p')o - no sin \{p + p')o] 
ifl'oC£;C[^'',^1ijO-Or(sin/c^ - sin Q sin\{k + /)o, 



Fig. 10 e 



[2 cos \ {k + Z)o cos + l)m cos |/o cos ^/Cm — COS COS |A;o] 

Fig. |T0^ = IglCit"-, t''}ij5mn[rs COS \{p + p')m " «7m sin |(p + p')m] 



X 



(A4) 
(A5) 
(A6) 
(A7) 



[2 cos + /)rra COS + On COS |/m COS — COS |/„ COS 

+ ifl'oC[^",^1ii^mnSini(/c + /)m[cB(Tmr(sinfcr - sinlr) + CECTmoisinko - sin/o)], 

where Si{pi) = sinpj; a is the gauge parameter (a = 1 yields Feynman gauge); go is the 
bare gauge couphng; the are (anti-Hermitian) generators of SU(3), normahzed so that 
trf^t^ = — (5"^/2; and the 7 matrix conventions are as in Ref. f^. One can verify easily that 
these rules reduce to the ones for the Sheikholeslami-Wohlert action when ( = 1, Ts = 1, 
and cb = ce- 



APPENDIX B: NUMERATORS 



It is convenient to sort the numerator -F(p, k) in Eq. (|4.9| ) according to the number of 
clover interactions: 

F = Fo + cbF,^ + ceF,^ + c|F,| + cbCeF^^^e + 4^4- (Bl) 

The Wilson-action term Fq is gauge dependent, but the other terms are independent of the 
gauge parameter. Below we give Fq in Feynman gauge [a = 1 in Eq. ( [A2D ]. The sum of 
internal quark and gluon momenta appears often, so let s = p + (p + /c) = 2p + k. From the 
Feynman rules, tedious algebra (performed by Mathematica and checked by hand) yields 

Fo{p,k) = 270 [sin(p + /i;)o{cosso - C^[rf - 1 + \{rl - l)s^]} + L{p + k) sinso 

+ <^7. [sm{P + kh (C'{coss. -[d-l + \{rl - l)s']} + 

i 

+ rsCL{p + k) sinsj 
+ sin So sin(p + k)o + r^C^ sin s ■ sin(p + k) 

- L{p + A;){cosso + eid - 1 - \{rl + l)s']}, (B2) 
where = I]j(2 sin |gj)^, and L{q) = 1 + moa + \rsC,q^ — cos go; 
FcB = i'TsC^^ili sin(p + k)j - 7j sin(p + k)i] sin ki sin \sj cos \kj 

— iC'^L(p + k)'^'yi sin ki cos ^kj cos + 51 ^« sin(p + k)i cos ^kj cos |sj; (B3) 
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Fcj, = ir.sC^'yo sin ko ^ cos sin ^Si sm{p + k)i 



— iC^7o sin |so cos |fco sin k ■ sin(p + fc) + sin koL^p + fc) ^ cos cos 

+ iC7 ■ sin k cos |A;o[sin(p + A;)o sin ^sq — + k) cos |so] 

— ir^C^ sin /cq sin(p + k)^ ^ 7^ cos |A;i sin |si 

i 

+ sin k ■ sin(p + k) cos |A;o cos |so + sin /eq sin(p + A;)o ^ cos ^fej cos |sj; (B4) 



= -iC^7o sin(p + k)o J2 sin^ h cos^ ifcj 

+ f C^7 ■ sin ^ sin /cj sin(p + k)i cos^ — |C^7 ' sin(p + A;) ^ sin^ ki cos^ |A;j 
+ ^ 7j cos^ ^ki sin fcj [sin(p + k)i sin /c^ — sin(p + k)j sin fcj] 

+ lC'^L{p + k)J2 sin^ fci cos^ ifc^-; (B5) 



cecb 



C 7o sin fco 2^ sin(p + sin fcj COS + f C sin(p + A;)o sin fco 2^ 7i sin fcj cos -^kj 



(B6) 



^4 = K [^70 sin(p + k)o - <7 ■ sin(p + fc) + L(p + A;)] x 



COS^ 1^0 X! ^^^^ ^« + ^^^^ ^0 X] '^^^^ 

(7 ■ sin k) [sin ■ sin(p + A;)] cos^ ^/cq + sin^ ko ^ 7j sin(p + k)i cos^ ^/cj 



(B7) 



APPENDIX C: INTEGRATION OVER LOOP MOMENTUM 



To apply Eqs. ( ^.351 ), ( ^^.201 ), and ( |2.26| ), one must obtain the self-energy functions Ap 
and C for real Pq and analytically continue them to imaginary values iE. The analytical 
continuation must be approached cautiously in the representation of the self-energy functions 
as integrals over the loop variable ko. The subtleties arise even with the continuum self 
energy. Setting the internal quark momentum to q, one has 



d'^q dm + i{d-2)4 



{2'kY {p- qf{q^ + m? 



(CI) 



In the complex go plane there are poles at ±iE{q) and at po ± i\p — q\- (Here, in the 
continuum, E{q) = y/q'^ + m^.) The integration contour is the real axis. If one sets po = 
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iE{p), however, pole at po — i\p — q\ crosses into the upper half-plane when |p — q^l < E{p), 
and the contour must be deformed to accommodate it. Numerical integration packages are 
not clever enough to deform the contour; it is usually the real axis by default. 

A better choice is to let q = p+k. In the complex fco plane the poles are at —po±iE{p+k) 
and ±z|fc|. Then setting po = iE{p) moves a pole across the real axis if |p| > |p + For 
the wave function and the rest mass, one wants external momentum p = 0, so the pole 
cannot, in fact, cross. When the quark mass is large, however, the pole comes very close to 
the real axis, producing a sharp peak near /cq = 0. Such integrals are difficult to estimate 
robustly. Even worse are the derivatives needed for the combination D, cf. Eq. (^.27^: 



partial differentiation with respect to pi and po each induces infrared divergences, which 
must, however, cancel in the total derivative. 

Consequently, it is prudent, if cumbersome, to integrate over ko analytically |^,|T2|,|13 



The easiest way is contour integration. In a lattice theory, one can proceed as follows, see 
Fig. |TT|a. One is to integrate on the real axis over the interval [— 7r,7r]. At ivr one extends 
the contour vertically to ±7r + ioo. Then one closes the contour with a segment from n + ioo 
to — TT + zoo, resulting in a rectangular contour. 

Since the integrand is a periodic function of ko, the contributions from the two vertical 
sides cancel. For the Wilson action, the contribution from the top vanishes; for improved 
actions, however, this need not be the case; then the top must be subtracted out explicitly. 
This can be done elegantly by the conformal mapping z = e^''°. The rectangular contour 



maps into the contour shown in Fig. pJ]b. The top of the rectangle maps into the small circle 
around the origin; whether it needs to be subtracted back out, or not, it is always correct 
to take only the outer circle at \z\ = 1. The final contour is shown in Fig. |ll|c. 

As a function of z the integrand has poles at physical locations and at the origin. The 
latter restore the contribution from a nonvanishing "top of the rectangle." It is then straight- 
forward to integrate over ko for all forms of ko dependence appearing in Eqs. (p2|)-| 



APPENDIX D: PAULI-VILLARS FUNCTIONS 

Here we give explicit expressions for the self-energy functions with the Pauli-Villars 
regulator. They are needed to reconstruct the lattice results from the coefficients in the 
tables. 

After introducing Feynman parameters and performing textbook manipulations one ob- 
tains 



sgUp'; rn) = dx (4m + 2zx^) In (Dl) 



where 

AA(a;) = xA^ + (l-x)m^-Fx(l-x)p^. (D2) 
Carrying out the integration over x and setting A — 0, one finds 



Ml 1 \ Cp \ 
J^'yi-m ;m = -7^ — 



1 - 2m^)(^(m^) + 1 - Inm , (D3) 
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FIG. 11. Contours for integrating over ko. (a) The complex ko plane; (b) conformal map of 
onto the complex z plane; (c) the final contour in the complex z plane. 
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and, recalling Eq. ( |4.14|) , 



where 



m 



m 



Q In + Vl — tanh ^ Vl — 
(^1 Inm^ — -y/4m2 — 1 tan~^ y/Am'^ — 1 



As m 0, (p{m ) ^ — 1 + Inm + m + Inm ); as m oo, (p(m ) 
tiating with respect to po, one finds 

loTT^ L 



where 



m ; m) 



167r2 



77(m2) + iln(AV"^^) 



^(m^) = m ^(1 — Inm^) — (1 — m ^)y9(m^) + ip{m'^), 



Tjijn^) = Inm^ — ipiju^) — ip^m^), 



with 



ipim^) = I Inm^ + 



(1 - 4m2)-V2 tanh"^ VT^^Am?, <\ 



{Arr? - 1)^^/2 tan"^ \/ Am^ - 1, > 



2^1 
4 



As 



0, ■ip{m'^) — > — m^(lnm^ + 1); as — > oo, ipiw?) \ \n.rn? . 
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